Ts = [12,15,20,23,26,30,35,40,50,60,80,100,130,160,190,220,250,275,296];
figure
hold on
dataall=zeros(19,1997);
for num=1:length(Ts)
    filename=['VV',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1),data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1);
    dataall(num,:)=data(:,2);
    ylim([20,500])
    title('VV')
end

figure
[X,Y] = meshgrid(Es,Ts);
pcolor(X,Y,dataall)

shading interp
colormap('jet')
clim([10,50])
xlim([505,525])
%%
figure
hold on
dataall=zeros(19,1997);
for num=1:length(Ts)
    filename=['VV',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1)/8.06554,data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1)/8.06554;
    dataall(num,:)=data(:,2);
    ylim([20,500])
    title('VV')
end
figure
imagesc(Es,Ts,dataall)
shading interp
colormap('jet')
clim([20,40])

figure


[X,Y]=meshgrid(Es,Ts)
mesh(X,Y,dataall)
colormap('jet')
clim([15,60])
xlim([55,70])
%% Imagesc plot
Tplot=0:305;
dataplot=zeros(length(Tplot),1997);
Tline=[];
for num=1:(length(Ts)-1)
    Tline(num)=(Ts(num)/2+Ts(num+1)/2)-0.1;
end
for idx=1:length(Tplot)
    numT=sum(Tplot(idx)>Tline)+1
    dataplot(idx,:)=dataall(numT,:);
end
figure(52);
imagesc(Es,Tplot,dataplot)
shading interp
colormap("jet")

clim([0,80])
xlim([500,530])
ylim([0,300])